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1-^ I We reply to the Comment of T. K. Kjeldsen and L. B. Madsen [e-print arXiv:physics/0508213] 

and acknowledge their criticism related to imperfect composition of our model 3ag molecular 
state [V. I. Usachenko and S.-I. Chu, Phys. Rev. A 71, 063410 (2005)]. However, we can- 
not agree with the authors suggesting this critique also as an irrefutable argument and/or 
evidence in support and justification of their opposite and incorrect (viz. inconsistent with rel- 
evant experiment) orientation dependence of N2 ionization rate calculated within the velocity 
^ ' go,uge version of molecular strong-field approximation. We disclose the actual reason of the 

?>^' contradiction and demonstrate that appropriately composed Sa^ state (modified according to 

the Comment's critique) rather confirms the alternative earler calculation [A. Jaron-Becker, A. 
Becker and F. H. M. Faisal, Phys. Rev. A 69, 023410 (2004)] suggesting correct orientation 
dependence of N2 ionization, contrary to respective results of the Comment's authors applying 

^ I the same approach and procedure of 3(jg composition. 
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1. INTRODUCTION. 

The current state-of-the-art in strong-field molecular ionization theory seems to be 
still far away from a sufficiently clear insight due to numerous controversial results ob- 
tained within different approaches and methods. As an appropriate illustration to mention 
here is the problem of strong-field ionization in laser-irradiated N2 molecule, which re- 
cently received a special consideration within both the velocity-gauge (VG) [1-5] and the 
length-gauge (LG) [3, 6] formulations of conventional strong-field approximation (SFA) [7]. 
In particular, the orientation behavior of N2 ionization rates was found in [3] to be quite 
different if calculated in the VG version versus the LG version of the molecular SFA ap- 
proach and/or molecular tunneling theory [8]. Based on the same self-consistent Hartree- 
Fock (HF) numerical procedure as previously proposed in [1, 2] for the 30"^ highest- occupied 
molecular orbital (HOMO) composition, the VG-SFA results [3] suggested a predominant 
A^2 ionization when the molecular axis is perpendicular to the incident laser field polariza- 
tion (see also Fig. 2 in the Comment [4]). Thus, in contrast to their LG-VGA results [3], 
the VG-SFA results of the Comment's authors suggest a counterintuitive behavior for N2 
ionization, which contradict to both relevant experiment [9] and alternative calculation 
[2]. Even though the Comment [4] applies in fact the same VG-SFA approach and numeri- 
cal (HF-based) procedure of 30"^ composition as earlier proposed in [1, 2], the Comment's 
authors carefully avoid of quoting the opposite results of Ref.[2] and/or clarifying the 
reason of the contradiction. They prefer instead to criticize our imperfect composition of 
model 3o"g molecular state [5] as inappropriately symmetrized and oversimplified. More- 
over, they even suggest their criticism (being physically correct and acceptable as such) 
as an indispensable explanatory reason of the contradiction. In order to demonstrate 
the invalidity of such an interpretation (as well as the related Comment's statements), 
we have to identify the actual reason of the contradiction by means of applying a more 
accurate molecular 30"^ state appropriately modified according to the Comment's critique. 



2. THE LCAO-MO METHOD AND ASSOCIATED MOLECULAR WAVEFUNCTIONS 
OF 3^3 STATE IN COORDINATE AND MOMENTUM SPACE. 
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According to the standard linear combination of atomic orbitals (LCAO) and molec- 
ular orbitals (MO) method, the molecular orbitals are the mathematical constructs used 
to model the multi-electron molecular valence shells (each of a fixed discrete binding en- 
ergy ^o = — J^"-* and respective number iVg^"^ of identical electrons, similar to atomic 
orbitals in atom). For 3(Tg HOMO of gerade and bonding symmetry in N2 molecule un- 
der consideration, the corresponding approximate two-centered single-electron molecular 
wavefunction can be composed, for example, as an appropriate superposition of scaled 
hydrogen-like 2pz and Is (or 2s) atomic orbitals: 
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which are centered on each of the atomic cores and thus separated by internuclear distance 
Ro- Here iV^^^'a) = 2, whereas 



S^ [Ro) = J dr<pf (r + Ro/2) ^f ^ (r - Ro/2) 



(3) 



denotes the respective atomic orbital overlap integral. 

Note that presently the "— " combination of 2pz states in Eq.(l) is used, which does 
provide the even spatial parity of the {2p) Sag HOMO. Thus, even though the "+" com- 
bination of scaled hydrogen-like 2pz atomic orbitals alone earlier used in the literature 
[5, 6] to reproduce Sag HOMO in N2 proved to be a surprisingly well working in ex- 
plaining relevant experiment [10, 11], we have to mention here the correct Comment's 
critical remark that such combination doesn't provide the even (gerade) spatial parity 
and bonding symmetry of the (2p) Sag MO. Moreover, according to the MO-SFA model 
[1], the oversimplified composition (1) would result in destructive intramolecular interfer- 
ence (see also [1-2, 5]) and accordingly high suppression in N2 ionization (contrary to the 
experiment [10, 11], which shows no suppression). To eliminate the mentioned deficiency, 
in our present consideration a more accurate composition of the Sag HOMO is applied 
taking into account some admixture of a comparable contribution from atomic s-states 
mostly required to provide a good agreement with experiment (see also Sec. 3 for details). 
Thus, in compliance with the second Comment's critical remark, the Sag HOMO is to be 
further approximated by the coherent superposition of few different MOs corresponding 
to separate contributions from atomic states of a specified orbital symmetry (viz., the 
scaled hydrogen-like Is, 2s and 2pz orbitals) under consideration: 
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with the weight coefficients Aj [j = 1,2,3) being the relative contributions {\Aj \ < 1) 
from respective atomic states and considered as variational parameters to be found from 
the equation for minimum of respective molecular binding energy, the value of which is 
put to be equal to the experimental value. 



The model single-electron molecular wavef unctions (l)-(2) allow for analytical rep- 
resentation of the Fourier transform (or respective molecular state corresponding to a 
definite value of final photoelectron momentum p in momentum space) for each of one- 
electron two-centered single molecular orbitals contributing to the Sag HOMO in N2 under 
consideration: 
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with the polar angle O-p of photoelectron emission relative to internuclear molecular axis 
(cos(^p) = (p ■ Ro) /pRo). Here 4"^ = -kI/2 = -I^""^ = - (zf Vaj)V2 is the bind- 
ing energy of n-th initial molecular discrete state under consideration (viz., /^^""s^ (A'^2) ~ 

15.58 eV) and Zj is the effective charge corresponding to " ejfective" long-range Coulomb 
model binding potential of respective residual molecular ion, while aj = jqq is j-th Bohr 
orbital radius of respective contributing atomic orbital. The analytical expressions for re- 
spective atomic orbital overlap integrals (3) are presented in [5], moreover, the Coulomb- 
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Volkov correction factor C (/t„) = i/lKnl^'^^ / E) " is also introduced to matrix elements 
(5)- (7) to incorporate the long-range Coulomb electron- molecular ion interaction in the 
final continuum states due to the VG version of SEA applied [1]. 

3. DISCUSSION. 

Fig.l demonstrates the angular dependence of molecular wavefunction (4) in mo- 
mentum space presented for different Sag compositions, either with or without taking 
a comparable contribution from atomic s-states into account. Let's note first that the 
angular dependence corresponding to contribution from 2pz states alone (Fig. la) is con- 
siderably prolate along the internuclear axis since the 2pz states are predominantly ionized 
along the internuclear axis owing to the presence of the geometrical factor cos 9-p in Eq. (5) 
(cf. also with Eq.(20) and respective Fig. la in [5]). This, particularly, illustrates that 
the spatial (gerade or ungerade) symmetry of respective {2p) Sag MO doesn't essentially 
affect the angular dependence (viz., nodal plane along the molecular axis) in momen- 
tum space. Thus, irrespectively of the sign in the combination of 2pz atomic states, the 
respective {2p) Sag MO suggests a predominant photoelectron emission along the inter- 
nuclear axis in momentum space. The sign in the combination mostly alters the absolute 
value of local maxima in respective angular dependence due to different {sine or cosine) 
trigonometric factor arising from intramolecular interference of contributions of ionization 
from two separate atomic centers [1]. In particular, the sin [(pat-Rq) /2] factor in Eq.(5) 
corresponding to the "— " sign results in destructive intramolecular interference and a 



suppressed molecular ionization (e.g. in O2 molecule relative to Xe atom, see also [1, 2, 
5], for details). 

Accordingly, if the currently used "— " combination of 2pz states is supposed to con- 
tribute to the Sag MO, then a coherent contribution from other (Is and/or 2s) atomic 
states should be also taken into equally detailed consideration as described by Eq.(4). The 
respective Figures Ic-ld (versus Fig. la) demonstrate that the contribution from s-states 
to the {ls2p) ScTg and (Is2s2p) Sag MO dominates in momentum space at low photoelec- 
tron momenta, in a qualitative accordance with stated in the Comment. Unfortunately, 
the Reply's restrictions don't allow us to provide the reader with an extensive and suf- 
ficiently detailed information about the respective molecular ionization rates, it will be 
given elsewhere in our further related publications. So far we only have to stress that 
the corresponding relative (weight) coefficients Aj at Is, 2s and 2pz states under con- 
sideration proved to have comparable values, although not necessarily of the same sign. 
Namely, the Sag compositions with Ais and A2p (as well as A^ and ^2^) of opposite signs 
only proved to provide no suppression observed in experiment [10, 11] for N2 ionization as 
compared to ionization of counterpart atomic Ar of nearly identical ionization potential 
ji-^r) r^ 15.75 e\/. In particular , the total ionization rates calculated for Ais and A2p of 
opposite signs in (ls2p) Sag and (Is2s2p) Sag MO proved to demonstrate no suppression 
(and even some enhancement) in N2 ionization relative to Ar well consistent with rele- 
vant experiment [10, 11]. Our Figs.lc-ld also suggest that {ls2p)Sag and {ls2s2p)Sag 
MO corresponding to the coefficients Ai^ and A2p of opposite signs are always noticeably 
prolate in momentum space along the internuclear axis, although to a considerably less 
extent as compared to the {2p)Sag composed of pure 2pz states alone. 

Anyhow, our Figs.lc-ld clearly illustrate that a more accurate consideration doesn't 
necessarily result in the Sag HOMO to predominate in momentum space along the di- 
rection perpendicular to the internuclear axis (as suggested, e.g., by Fig.4e of the first 
reference in [3] and/or Fig. If in the Comment [4]). Since the generalized Bessel function 
are commonly maximized when photoelectron momentum p is parallel to the laser field 
polarization e, the ionization from {ls2p) Sag and {Is2s2p) Sag presented in our Fig.l is 
expected to be predominant when the A^2 molecular axis is aligned along the laser field 
(Ro||e). This difference thus invalidates the Comment's statement that incorrect (viz., 
inconsistent with experiment) orientation dependence of A^2 ionization rate the authors 
also found in [3] is just explained by the VG version of molecular SFA they applied. 
Our Fig. 2 well illustrates the above conclusion and displays the orientation dependence 
of total N2 ionization rates (see Eq.(13) in [5]) corresponding to different Sag composi- 
tions. As expected, the compositions of Sag HOMO taking a comparable contribution 
from p-states into account always suggest a predominant ionization of N2 for parallel 
orientation of internuclear axis relative to the laser field polarization. Whereas, the com- 
position of Sag accounting for a contribution from pure s-states alone is well seen to only 
give a predominant ionization of A''2 for perpendicular orientation of internuclear axis as 
suggested in the Comment. Such a behavior may thus arise only due to the contribut- 
ing Is or 2s states, which are spherically symmetric and equally well ionized along any 
spatial direction. Therefore, the form of respective nodal plane relative to the molecu- 
lar axis suggested by Eqs.(6) and (7) is solely determined by the cosine trigonometric 



factor cos [(pat-Rq) /2], which alone suggests the angular dependence to be oblate along 
the internuclear axis in momentum space (Fig. lb). The resulting orientation dependence 
with maximum ionization for the perpendicular orientation of molecular axis is thus quite 
naturally expected for ionization of H2 and H2 (with a single (Is) lag MO composed of 
pure s-states). For H2, such orientation dependence was recently revealed [12] using a 
tunnelling model with a Heitler-London type valence bond wavefunction. 

To summarize, our argumentation based on a more accurate consideration of Sag 
HOMO (4) is also qualitatively well consistent with alternative calculation [2] using the 
same VG-SFA approach as applied by the Comment's authors [3, 4]. Meantime, we also 
have to stress that opposite orientation dependence of N2 ionization presented in Fig. 2 
of the Comment [4] is based on the same numerical HF-based procedure of ScXg HOMO 
composition as applied in [2] under precisely the same situation (cf. with respective 
Fig. 7b in [2]). Therefore, we believe that orientation dependence of N2 ionization rate 
presented in our Fig. 2 is correct and its consistence with relevant experiment [9] not 
just an accidental. By means of direct comparison of our Fig.lc (and/or Fig. Id) with 
respective Fig.4e presented in [3], one can identify that opposite behavior of A^2 ionization 
suggested in [4] (as well as earlier in [3]) stems from a different angular dependence of 
3(7g molecular state in momentum space. Our Fig.l also suggests that such difference 
may only arise owing to relative contribution from s-states, which seems to be somehow 
overestimated (compared to p-states) within the particular (12s7p)/[6s4p] basis chosen in 
the Comment for numerical HF-based procedure of 3ag composition in coordinate space. 
This, particularly, becomes especially evident from Fig. 3 displaying the correct orbital 
shape for the Sag wavefunction in coordinate space presented in Ref.[13] (see Fig. 12 
therein), which obviously differs by a remarkably prominent domination of p-states from 
the resulting 30"^ state of substantially different shape and detailed structure suggested 
in [4] (cf. respective Fig.lc therein). Meantime, the numerical procedure applied in [3, 
4] for 3ag composition doesn't seem to be sufficiently transparent with respect to the 
extent of such relative contribution, which is thus hardly amenable to reliable analysis 
and/or available for direct comparison [14]. Note also that a composition overestimating 
the relative contribution from s-states in coordinate space is not suitable to Sag state in 
N2, which is known to have a different and more complex structure well dominated by 
p-states [15, 16]. 

Finally, we also have to clarify the main message of a different (although, somewhat 
related) research [17] the Comment's authors quite inappropriately and very ambiguously 
invoked in support of their related statement suggesting the LG version to be superior to 
the VG version of the SFA under discussion. Recall that the cited similar conclusion of 
Ref.[17] (see also Fig.l therein to make certain) applies solely to ionization of atomic states 
of only odd spatial parity, such as p-states, which are not symmetric with respect to the 
inversion of all electron coordinates to opposite direction. Whereas, for the molecular Sag 
state (which is known as surely gerade, i.e. of even spatial parity), the VG and LG versions 
of SFA, according to [17], are expected to give the same results, contrary to respective 
VG-SFA findings of Ref.[3] and the context of main related statement of the Comment [4]. 

4. CONCLUDING REMARKS. 



The criticism of the Comment [4] related to our previous imperfect composition of 
the Sag molecular wavefunction [5] was assessed to be physically correct as such, however 
not affecting essentially our related final results and conclusions based on the velocity 
gauge (VG) version of molecular strong-field approximation (SFA). In particular, when 
applying an appropriately modified Sag wavefunction, we demonstrate that orientation 
dependence of N2 ionization is to be correct (viz., consistent with experiment) and thus 
the Comment's critique cannot be accepted as an irrefutable argument and indispensable 
evidence in support of the authors' opposite respective VG-SFA findings [3, 4]. This, 
particularly, invalidates the Comment's main statement about unavailability of the SFA 
in its VG formulation for adequate description of molecular strong-field ionization when 
an accurate molecular orbital is applied. Our conclusion is also well consistent with re- 
cent Erratum [18] suggesting in fact analogous corrections to analytical expressions for 
molecular ionization rate previously numerically calculated in [1, 2] for diatomic molecule 
having a HOMO of gerade symmetry. Those corrections were also confirmed in [18] as not 
affecting their correct final results and conclusions discussed elsewhere earlier including 
the interference suppression of molecular ionization and orientation dependence of N2 
ionization [1, 2], which proved to be quite opposite to respective VG-SFA result of the 
Comment applying the same approach and pure numerical HF-based procedure of Sag 
composition. 
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FIGURE CAPTIONS 



2 

of Sag HOMO wavefunction in 



1. Fig. 1 (Color online). The squared module $3o-g(pAr,RQ 
momentum space (the Fourier transform) plotted vs the angle 6p of photoelectron emis- 
sion (relative to the internuclear molecular axis lined up in the horizontal direction) and 
calculated for four different Scg compositions: (a) (Is) Sag MO composed from Is states 
alone, (b) (2p) Sag MO composed from 2pz states alone, (c) (ls2p) Sag MO composed 
from Is and 2pz states, (d) (Is2s2p) Sag MO composed from Is, 2s and 2pz states. These 
angular dependencies are all presented for photoelectrons only emitted along the polar- 
ization of Ti:sapphire incident laser field of the wavelength A = 800 nm {fiuo = 1.55 eV) 
and fixed intensity / = 2 ■ IQ^^W/cm? due to absorption of different number A^ of laser 
photons beginning from minimum one A'"o = 18 required for ionization of A^2- 

2. Fig. 2 (Color online). The total N2 ionization rates corresponding to Sag compositions 
presented in Fig.l and plotted vs the angle G of the internuclear axis orientation relative 
to the polarization of Ti:sapphire incident laser field of the wavelength A = 800 nm 
{huj = 1.55 eV) and fixed intensity / = 2 ■ lO^^W/cm'^. The ionization rates are all 
normalized to respective ones at the angle = and presented for four different Sag 
compositions: {ls)Sag MO composed from Is states alone (solid line); (2^)30-3 MO 
composed from 2pz states alone (dashed-dotted line); {ls2p)Sag MO composed from 
Is and 2pz states (short dashed line), {ls2s2p)Sag MO composed from Is, 2s and 2pz 
states (long dashed line). This figure is to be compared with Fig. 2a of the Comment and 
Fig. 2a of relevant experiment [9]. 

3. Fig. 3 (Color online). The orbital shape of the correct Sag HOMO in N2 molecule im- 
aged by recording the high harmonic spectra and presented in [13] (see Fig. 12 therein 
and related explanation). Note that this is not the orbital density, but a wavefunction 
to within a global phase. This figure is to be compared with respective Fig.lc of the 
Comment displaying the alternative Sag HOMO numerically composed within particular 
(12s7p)/[6s4p] basis set. 
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